Strong Soret effect in one dimension 
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£NJ . We consider a one-dimensional gas of two kinds of particles with different masses interacting 

» I ' through short range interactions. The system exhibits an extreme form of the Soret effect: when 

f> |_ the ends of the system are in contact with thermal baths of different temperatures, there is complete 

■ separation of the species. We show how this separation can be well described in the Boltzmann 

approximation and discuss the origin of this odd behavior. 
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i I. INTRODUCTION 

. When a binary fluid is subjected to a temperature gradient, the densities p + and p_ of the species become space- 
dependent in such a way that the relative concentration {£)/ p- [x) changes along the direction of heat flow. This 
I \ is the Ludwig-Soret, or Soret effectta. That this effect happens does not in itself need an explanation: once the 
temperature difference takes the system away from equilibrium and breaks translational symmetry, there is no reason 
why the ratio of concentrations should stay constant. _ 

However, the Soret effect has practical applications in the separation of speciescl, and this requires apquantitative 
prediction of its magnitude in-,each case. The*e have been several analytic approaches to do this for gasestil, condensed 
phasesa, grains in suspension^, porous mediaH, etc. (For a review, seeQ). 

From a purely theoretical point of view the Soret effect is interesting because it cannot in general be discussed in 
terms of local equilibrium, in which a space-dependent local temperature fully accounts for the concentration ratio 
- as for example the local pressure accounts for the variation of the density in a gas column under gravity. To 
see this explicitly, we consider here a model belonging to a family of systems for which such a local equilibrium 
approximation yields strictly zero effect, contrary to observation. They consist of two species of different masses m+ 
and m_, interacting through a pair potential V(ri — Tj) which does not depend on the particle type (the masses of 
t-H particles i and j). A short calculation shows that the partition function Z of such a system is of the form: 
> ; 

— lnZ(p+,p_,T) = -p + \n P+ - p_ lnp_ + -p + ]n(m+T) + -/?_ lu(m_T) - (3F(p,T) (1) 

cn : v 

where F depends on p + , p_ only through their sum p = p + + p_ , and d is the dimension. Assume now that under 
heat flow we can divide space in small cells 5 d (x) with densities p+(x) and p^(x), within which the system is in 
equilibrium at temperature T(x). The total partition function will then be a product expressions (|) for each cell. If 
we now maximize it with respect to p+{x) and p~(x) (with fixed number of particles) we get: 
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g; a-=-^ — (p,T(x))-ln(p_) + -ln(m-) (2) 

with a± space-independent Lagrange multipliers. Subtracting these equations we get p + / p- = constant, i.e. no Soret 
^ 1 effect. 

In this paper we study a very simple one-dimensional model which belongs to the class described above, and has a 
most extreme form of Soret effect: there is total phase separation even for arbitrarily small temperature differences. 
The two pure phases are separated by an interface in which species are mixed, its width is a function of the parameters. 
Although this full separation is not realistic, it is interesting to exhibit a model which, for large sizes, is driven far 
from equilibrium by an arbitrarily small difference of temperature. 

We also study in this model the closely related question of the laws for heat transfer. In regions in which phases are 
pure this transfer has an anomalous spatial and temperature dependence, just as found in oscillator chainsu. Within 
the interface we show that the behavior becomes the usual Fourier law, as recently found by Dhara in a model similar 
to the present one, but in which particles are forced to stay mixed. 
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II. THE MODEL 



We shall consider a one dimensional system of length L consisting of iV point particles of two species: N + heavy 
and N- light particles of masses m+ and m_, respectively. 

Particles interact through an infinitely narrow potential of height V, so that when two particles meet they collide 
if their center of mass energy 

1 m 1 m 2 , 2 . . 
E c =- ■ [V1-V2) <V 3 

2 mi + m 2 

and ignore each other otherwise. The collision between particles conserves energy and momentum: 

vi — ► Ui = vi H (v 2 - vi) (4) 

mi + to 2 

"2 >V 2 = V 2 ; («2-«l) (5) 

mi + m 2 

These collision rules are reversible and satisfy detailed balance. Between collisions the particles evolve freely across 
the system. Particles colliding against the walls x = and x = L rebounce with a velocity with random distribution 
corresponding to thermalisation at temperatures T c and X^, respectively: 

m / mv 2 \ 

P a (v)dv = -^uexp \ --^r J dv a = c,h (6) 

This guarantees that if T c = T), the system equilibrates to the Maxwell distribution. 

In fact, only the potential between particles of different species is relevant, since particles of equal mass only exchange 
their velocities in a collision, which then just amounts to exchanging their labels — a one dimensional peculiarity. As 
a consequence of this, we can consider this model as having a potential V(ri — Tj) which is independent of particle type 
(and hence belongs to the family (|l|)), and at the same time consider that particles of the same type are transparent 
to one another, and hence our model is a particular case of the one introduced by Widom and RowlinsorEB. 



III. PHASE SEPARATION 



We use a simple molecular dynamics in the simulation, calculating the minimum collision time and letting the 
system evolve freely between collisions. We perform simulations for different number of densities, and different values 
of the length ranging from L = 100 to L — 2000, which allows us to consider finite-size effects. 

We have checked that in the equilibrium case (the temperature in both sides of the system are identical) the two 
species are homogeneously distributed in space and their velocity distribution is Maxwellian. We also have checked 
that for the case of one type of particles with two different temperatures, the average energy is the one corresponding 
to an homogeneous temperature \JT c Th- This is because if we consider only one type of particles, in the elastic case 
the particles simply cross the system without any interaction with other particles and it is easy to verify that the 
distribution of velocities is just composed of two Maxwellian distributions at different temperatures for particles flying 
in each direction. 

In Fig|l] we show the evolution of the particles after a long thermalisation time, starting from a configuration with 
random positions [L = 100, V/T c = 10, T = Th/T c = 10, m_/m + = 0.2). We can see that the space is divided in 
two regions with only a single species, separated by an interface in which both species mix. Away from the interface, 
because all particles have the same mass, the system can be considered as non-interacting. Hence we have the following 
picture: the heavy particles collide with the hot wall at x = L. They fly unperturbed (modulo relabelings) throughout 
the the heavy-particle phase, until they meet the first light particles around the interface. After one or more crossings 
and collisions with light particles, they return to the heavy-particle phase, where they are free again. The same can 
be said of the light particles on the left. 
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FIG. 1. Plot of the time evolution of the positions of the particles. The large circles correspond to a particles with m+ = 1 
while the small ones to a m_ = 0.2. The walls are a T c (x = 0) = 1 and Th(x — 100) = 10. The interaction potential is V = 10. 
There are 20 particles of each type. 

A simple computation helps clarify what happens in the interface. Consider a system with a light and a heavy 
particle. Estimating their velocities according to the temperatures of the walls, we have: 



If the light particle is on the left and the heavy one on the right: 



(7) 



which gives a typical center of mass energy: 



IE. 



(1 ) _ (T c + T h )(m+ + m_) - (^T c m- + ^/T h m + f 



CM 



(m + + m_) 



(8) 



• If the heavy particle is on the left and the light one on the right: 



(9) 



giving the typical precollision center of mass energy: 



(2) _ (T c + T h )(m + + m_) - (y/T h m- + ^T c m + ) 2 



2E cm - 



(m + + m_) 



(10) 



In fact is easy to show that: 
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(i) 

CM 



2E. 



(2) 
CM 



< 1 if m + > to_; T c < Th 



(11) 



Hence, we see that the probabilities of collision are not symmetric: while this does not prove that phase-separation 
exists, it gives us the sign of the effect. Note that there is a similarity with asymmetric exclusion processes^, although 
unlike that case here detailed balance holds in the bulk. 
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IV. INTERFACE 



In Fig|2| we show the time-averaged particle densities p+ (x) and p_ (x) : 

p±(x) = 1 It \im [ dt'y^S{x-Xi{t')) (12) 

for a system with the same parameters as Fig|l|. We see, as before, that the two types of particles are completely 
separated, but the interface seems much thicker than in Fig.|l|. The reason for this discrepancy is simple: the position 
of the interface fluctuates with time. Upon time-averaging we are actually measuring the instantaneous interface 
width convoluted with the dispersion in its position. 
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FIG. 2. Plot of the densities p+ and p~. The parameters of the model are the same as in Fig.|l| Data obtained with 100 bins 
taking measures every 10 collisions during a time corresponding to 500.000 collisions per particle. The diamonds correspond 
to the heavy particles while the circles correspond to the light ones. 

The fluctuations in the position of the interface AL scale with the system size L. In Fig.|| we check the assumption 
that, just like in an ordinary equilibrium system AL oc \f~L~. In order to do this, we plot the product of densities 
p+(x)p-(x), a quantity that is large only in the interface, for several system sizes and we verify that, indeed, the 
function scales with \[~L. In || we have considered again the same parameters as in Figs.[l] and ||, and systems lengths 
from L = 100 to L = 800. Rescaling the averaged interface with y/~L, the curves collapse. 
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FIG. 3. Plot of the product p+(x)p-(x). For different sizes L, centering and rescaling the interface with yL. The curves 
correspond to L = 100, 200, 400, 800. 

In Figures ^ we plot p+p- for a system of length 100 and V — 10. On the left figure we vary the temperature of 
the hot wall keeping the mass ratio constant, while on the right figure we have changed the mass ratio while keeping 
the temperatures constant. These plots show us the shift in terms of the parameters of the average location of the 
wall. 
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FIG. 4. Product p+(x)p-(x) for different values of the parameters. Left: Th/T c = 4,5,8, 10, 15 and m_/m+ = 0.2. Right: 
m„/m+ = 0.05, 0.1, 0.15, 0.2, 0.25 and T h /T c = 10. 

As we seen above, the time-averaged width of the interface is dominated by the time fluctuations of the wall position, 
rather than by its instantaneous width. Hence, it is convenient to define a parameter that reflects the mixing of the 
particles at every time, thus giving a measure of the instantaneous width of the interface. For each heavy particle we 
count the number q t (= 0,1, 2) of its nearest neighbors that are light particles. Our parameter is then: 
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(13) 

For a sharp interface 9 = 0, and for a very mixed system q is of order N. In FigJE] we show the frequency distribution 
of q averaged over many time-steps for various sizes of the system. We see that the averaged distribution of q (which 
reflects the instantaneous width of the interface) is independent of the system size, as expected. We can hence study 
how the width changes with the various parameters, in a system-size independent way. 
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FIG. 5. Check of the size independence of the averaged frequency of q. The parameters are V/T c — 5, Th/T c — 10, 
m_/m+ = 0.2. 



On (a) we can see how the sharpness of the interface increases as we decrease the mass ratio. The same effect 
occurs when we increase the interaction potential (b) or the temperature ratio (c). The physical interpretation of 
these tendencies is easy to understand. The system is completely mixed when the two thermal walls are at the 
same temperature. In fact, in this case the asymmetry ( |ll| ) in the probability of collisions disappears (A = 1), and 
equilibrium is reached. Obviously if the masses are equal or the interaction potential is zero, the system is also 
completely mixed because there are no collisions. 

We have found no evidence of a threshold in the values below which there is no separation for an infinite system. 
Indeed, the limit of 'mixed particles' (e.g. at equal temperatures or equal masses) is only achieved when the interface 
length becomes larger than system size. 
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FIG. 6. Dependence of the distribution of q on the parameters, (a) Th/Tc — 10, m_/m+ = 0.2 and V/T c = 4,5,6,7; (b) 
Th/Tc = 10, V/T c = 5 and m_/m+ = 0.05, 0.1, 0.2, 0.3 (c) m_/m+ = 0.2, V/T c = 5 and Th/Tc = 5, 10, 15, 20 



V. BOLTZMANN EQUATION 



Let us now write down a Boltzmann equation for the probability distribution function of the velocities of the 
particles, and check that it reproduces correctly the phase-separation effect. 
We have to deal with two coupled equations: 



df{x,v,t) 

at 



M 

dx 



\v - u\(f(x,v',t)g(x,u',t) - f(x,v,t)g(x,u,t))du 



(14) 



dg{x,v,t) 
dt 



dg 
dx 



\v - u\(f(x,u',t)g(x,v',t) - f(x,u,t)g(x,v,t))du 



(15) 



where f(x,v,t) and g{x,v,t) are the probability distributions for the two types of particles. As usual, u' and v' 
are velocities the particles should have before the collision in order that the velocities after collision are u and v 
respectively. These equations must be solved with the adequate boundary conditions @. In this equations we have 
omitted the interaction potential, in fact there would be a O-function which guarantees the collision. 

The approximation that is made in equations ( |l4| ) and ( |i"5| ) is that the correlations between particles are neglected: 
each particle collision is independent of the collisions the particle has suffered before. More formally, this means that: 



f 2 (x, x, v, v', t) = f(x, v, t)f(x', v', t) 
g 2 (x, x', v, v', t) = g(x, v, t)g(x', v' , t) 



(16) 
(17) 



Unfortunately, as far as we know, the system ([uj) and-/[l5|) cannot be solved analytically. We solve it numerically 
using a stochastic process based on the Bird algorithms which has been provedt^l to converge to the solution of the 
Boltzmann equation. We have used this algorithm for this system taking into account the delta-like potential. 

In Fig.^ we show the profile of the spatial density of the two types of particles for a system consisting of 200 particles 
of each type for a potential V/T c = 4, mass ratio m+/m_ = 0.2 and temperature ratio Th/T c = 10. We can see that 
the two normalized profiles are almost identical, which confirms the validity of the Boltzmann approximation. This 
good agreement holds for any value of the parameters. 
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FIG. 7. Left: Densities profiles (left) and p+(x)p-(x) (right) for a system of length 100 with 200 particles of each type 
and V/T c — 4, m_/m+ = 0.2 and Th/T c — 10. The squares represent molecular dynamics while the circles the results of the 
Boltzmann equation. 

The fact that the Boltzmann equation is a good approximation, at least away from the interface, is easy to un- 
derstand. As the two species are separated, the most likely situation after a collision is that at least one of the two 
colliding particles goes freely to its corresponding wall, it loses all memory of the correlation with other particles. 
Thus, the neglect of the correlations implicit in the Boltzmann equation is a good approximation. 

A special note should be made about the interface fluctuations. The Boltzmann equation has a unique solution for 
every finite length L (it is first order). However, if we were to compute the stability matrix around the solution, we 
would find that the mode corresponding to translation of the wall becomes softer as L — > oo. This in turn means that 
corrections to the Boltzmann approximation become more and more important in that limit, and they yield the wall's 
fluctuations. (This question becomes more familiar if we bear in mind the analogy with a 2D Ising ferromagnet with 
the spins on right boundary set to +1 and on the left to —1. The Boltzmann equation is analogous to a mean-field 
solution, which places the interface in the middle, and neglects the fluctuating displacements of the domain wall.) 



VI. HEAT TRANSPORT - RECOVERY OF FOURIER LAW FOR BROAD INTERFACE 

In one dimensional systems, the obtention of Fourier's law is still an open questionS As mentioned above, this 
one dimensional model has the peculiarity that particles can be considered non-interacting everywhere except in the 
interface, where both species mix. This reflects itself in an extremely pathological form of heat transfer: the only 
spatial gradient in kinetic energy of the particles occurs across the interface, since far from the interface we have a gas 
of noninteracting particles. It thus seems reasonable, if we wish to recover the usual laws for heat transfer, to place 
ourself in conditions such that the interface is broader than the sample itself and species are mixed. 

In order to study the temperature profile in the interface we have chosen very small temperature differences and 
species with very similar masses. In Fig. 8 we show the energy profile for a system of length L = 100 with 100 particles 
of each type, with masses m+ = 1 and m_ = 0.95, respectively. We plot the energy profile for m + for different 
values of the temperature ratio, Th/T c = 1.05, 1.10, 1.15, 1.20. We can clearly see that the energy distribution changes 
linearly along the system, and that it is proportional to the thermal gradient. In fact, the straight lines fitted have 
all the same slope, indicating that the conductivity in the interface is the same for all cases. 
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FIG. 8. Plot of the energy profile in function of the temperature ratio along the system for the particles with mass m+ = 1 
when the other particles have a mass of m_ = 0.95. The length of the system is L = 100. The temperature ratios plotted are, 
from top to bottom , Th/T c = 1.20, 1.15, 1.10, 1.05. All the fitted lines have the same slope. 

Even in this small gradient case we can see that the mass distribution is not uniform. In fact, the massive particles 
are located closer to the hottest wall while the lightest ones are closer to the coldest wall. This fact can be seen in 
Fig. 9, where we have plotted the normalized particle densities for the same case as in Fig. 8 for Th/T c = 1.30. Note 
that the mass profile is also linear in small thermal gradients. 
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FIG. 9. Plot of the mass profile along the system for the particles with mass m+ = 1 when the other particles have a mass of 
m_ = 0.95. The length of the system is L — 100. The temperature ratio considered is Th/Tc= 1.30. The curve with positive 
slope correspond to m+ = 1. 



The fact that we have found the Fourier Law for the broad-interface case is not surprising, in view of the similar 
result of DharQ mentioned in the introduction. Note, however, that unlike the present situation, in Ref.EI particles 
with different masses are never allowed to cross, so the Soret effect is avoided completely. 
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It can also be seen that for a fixed small thermal gradient, the temperature profile depends linearly on the difference 
of masses. This fact suggests the possibility of expanding the Boltzmann equation in powers of this parameter and 
making an estimation of the conductivity for the case in which the masses of the two species are very similarE-3 . 



VII. CONCLUSIONS 



In this paper we have numerically solved a model with two types of particles with different masses. We have shown 
that a thermal gradient is enough to separate completely the two species, at least if the system is large enough to 
accommodate the interface. The phase separation is due to the collision asymmetry ( |ll| ) provoked by the thermal 
gradient applied. For large enough systems, the heavy particles are closer to the hottest wall while the light particles 
are closer to the coldest one. Between these two regions there is an interface where the two types of particles live 
together. In the regions occupied by only one type of particles the kinetic energy gradient is zero because the system 
is a noninteracting gas, so that all the energy density drop is localized in the interface (a finite region of the space). 
We have also studied the behavior of the interface for different values of the characteristic parameters of the system, 
showing that its instantaneous width is independent of the system size. By means of a numerical solution of the 
Boltzmann equation for this systems, we have checked that the total phase separation is also obtained within this 
approximation (equations (fil|) , (|l5|) ) . We have also investigated the limit of 'broad interface', in which the interface 
is broader than the system itself, and shown that Fourier's Law is recovered in that case. 

To sum up, we have introduced a very simple one dimensional model which has an extreme Soret effect for an 
infinite system, while in a finite system reproduces not only Fourier's Law but also a linear distribution of densities 
along the system. 
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